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Preface 

This final technical report consists of three documents organized into three 
chapters. 

The first chapter is entitled "A Semi-Analytic Monte Carlo Radiative Transfer 
Model for Oceanographic Lidar Systems". This chapter summarizes the computer 
model SALMON as it was developed for homogeneous media. It discusses the advantages 
of the model over more conventional Monte Carlo techniques used by other researchers. 
Chapter 1 served as a preliminary version for a paper published by L. R. Poole, 

D. D. Venable and J. W. Campbell. 

The second chapter is entitled "Sensitivity of Airborne Fluorosensor 
Measurements to Linear Vertical Gradients in Chlorophyll Concentrations". This 
chapter gives detailed discussions on the application of the stratified model 
to useful measurement related problem. Much of the mt. ..al introduced in chapter 
2 is expanded in chapter 3. This chapter serves as a preliminary paper for 
proposed publication. 

The final chapter is entitled "SALMON Modifications for Non-homogeneous 
Media". This chapter discusses the extension utilized to make the SALMON code 
appliable to stratified media. Some comparisons of the stratified and homogeneous 
models are discussed. 

In addition to the work presented here, several papers have been published 
by L. R. Poole et. al. concerning SALMON. The results contained in these reports 
were developed collabcrativeiy by D. D. Venable, L. R. Poole, J. W. Campbell and 
A. R. Punjabi. 
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Chapter i 


A Semiana lytic Monte Carlo Radiative Transfer 
Model for Oceanographic Lidar Systems 


L. R. Poole, D. D. Venable and J. W. Campbell 




Much attention has been focused recently on the application of Ildar 
techniques to oceanographic research problems (ref. 1), particularly those 
of Interest to biological oceanographers. Variations In time and space of 
Important biological parameters present a formidable obstacle for the 
researcher using conventional shipboard sampling and measurement techniques. 
Measurements made using shipboard or towed Ildar systems, such as the proposed 
Submerged UDAR and SONAR (SLAS) system, would* enhance conventional measure- 
ments, but would still be faced with the problems associated with spatial and 
temporal variability. However, It Is felt that carefully planned Integration 
of precise, but localized shipboard measurements with less precise, but synoptic 
measurements made using Ildar remote- sen sing techniques from aircraft or heli- 
copters could lead to a much greater understanding of Important research areas 
such as, for example, phytoplankton ecology (ref. 2). In addition, such a 
combination of measurements could serve as a surface-truth benchmark calibration 
of any future satellite-borne remote sensors addressing the same problems. 

Presently, oceanographic Ildar measurements are, at best, semi-empirical 
In that only a rudimentary understanding exists of'the effects of a number of 
factors (l.e.: attenuation of the laser beam In the water medium -diffuse or 
beam attenuation, or vertical Inhomogeneities In the medium). A technique Is 
needed for Investigating the accuracy of LIDAR system measurements over a wide 
range of environmental conditions, as well as for various boundary conditions 
and source/sensor geometric configurations. This letter presents a semi-analytic 
Monte Carlo radiative transfer simulation model (SALMON) which Is particularly 
well suited for addressing these fundamental questions. Sample results are 
presented which highlight the advantage of using SALMON as opposed to purely 
stochastic Monte Carlo techniques. 
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Monts Carlo simulation techniques have been used rather extensively In 
recent years to study radiative transfer In the atmosphere and oceans. In Its 
simplest digital form* a Monte Carlo model Is a purely stochastic construction 
of an ensemble of photon trajectories through the medium of Interest, with the 
length and direction of each trajectory segment governed by probability density 
functions derived frcm the basic scattering and absorption properties of the 
medium. In most applications, some technique or combination of techniques for 
reducing the variance In the computed results Is used In order to obtain meaning- 
ful results at a reasonable computer cost. Examples of Monte Carlo studies 

i 

of oceanic radiative transfer can be found In Plass and Kattawar (ref. 3), 

Funk (ref. 4), Gordon et al. (ref. 5), and Thomas and Guenther (ref. 6). The 
last of these examines, for the case of an unlimited fleld-of-vlew receiver, 
the Influence of scattering and absorption mechanisms on the performance of a 
bathymetric Ildar system, one which Infers water depth in coastal regions by 
transmitting a laser pulse Into the water and monitoring the return pulse from 
the water-bottom Interface. Thomas and Guenther state that a topic of major 
Interest not included In their results Is the quantitative effect of a limited 
receiver field of view. 

An obvious technique for limiting the receiver field of view In a Monte 
Carlo simulation of any Ildar remote-sensing system Is to accept as "successful" 
only those photon trajectories which terminate at a point lying within a pre- 
defined receiver aperture (at some height above the air-water Interface) and 
whose final segment Is In a direction lying within the confines of a specified 
angular field of view (see fig. 1). For most physical systems, the true 
geometric constraints are very restrictive, and, as a result, the probability 
of "success" for any given photon trajectory Is very small. The only choices 
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available to the researcher using a purely stochastic approach are to 
construct a very large ensemble of trajectories, thereby consuming considerable 
computer resources, or to relax the geometric constraints to more practical 
values and compute an average value of the detected signal over a mjch broader 
field of view than that associated with the actual sensor. 

A semi -analytic approach, however, one In which the detector signal Is 
evaluated through a combination of stochastic and analytic techniques, offers 
an effective means for modeling the Ildar system with realistic geometric 
constraints. The approach used In developing SALMON Is frequently called 
the "method of expected values" or the "method of statistical estimation'! 

(ref. 7). Consider again figure 1; suppose that the j— segnent of a photon 
trajectory constructed using standard Monte Carlo methodology Is In direction 
£ and terminates at some point q within the "Interaction volume", that 
voluma enclosed by an extension Into the water (with refraction effects) of 
the cone defining the tetector fleld-of-vlew. 

At this point In Its history, we can conceive of the photon as repre- 
senting a large number of photons, each traveling In the same direction and 
having the same amount of energy (the sum of which equals the energy level of 
the true photon at point “q). VI th a knowledge of the probability density 
functions governing photon motion, we can estimate the fraction of this large 
number of photons which, after Interacting with, say, a sediment particle at 
point q, (1) will scatter through an angle 0" toward a direction £' which 
lies within the small solid angle Aft subtended by the detector aperture 
A # p Sit point q, and (2) will proceed from point “q 1 to the aperture without 
further Interaction with the medium (or with the atmosphere, should we have 
to consider Its effects, also). Neglecting atmospheric effects, the estimate, 
or expected value, of the fraction of photons collected by the remote receiver 
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upon a scattering event at point q can be written as 

E ■ J> - iST ) " . 4n, * ed T $ o) 

where pie') Is the scattering phase function (assumed constant ever the 
small solid angle Aft), which gives the rate at which energy Is scattered 
Into an element of solid angle Aft about direction 8" (ref. 8). The term 

*cd 

• Is the probability that the photons scattered through angle d' will 

then transmit from point q to the air-water Interface (a distance d') with 

no further Interactions. The term c Is of course the beam attenuation 

coefficient for the water medium and T Is the Fresnel transmittance of 

s 

the air water Interface. An analogous relationship can be written for 
fluorescent emission from a chlorophyll-bearing particulate, with the scattering 
phase function replaced by a function characterizing the angular emission 
pattern of the particulate (usually assumed to be Isotropic). Upon each 
scattering event (or fluorescent emission) occurring within the interaction 
volume, the expected value of the energy reaching the detector Is recorded, 
and the total energy of the fraction remaining at point ~q Is reduced by the 
same expected value. The length and direction of the j+1 st trajectory 
segment are then selected using standard probability density functions, except In 
the case In which the. new segment Is In a direction toward the air-water 
Interface. In such Instances, the segment length Is drawn from a 
distribution which is truncated (ref. 9) to ensure that the photon will remain 
within the water medium and potentially contribute again to the cumulative 
detector signal. SALMON also makes use of the statistical weighting technique 
(ref. 9) whereby photons are prevented from being absorbed by the water 
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molecules themselves (and thus removed from the field), but rather allowed 
to proceed to further Interactions, with an appropriate adjustment In 
energy level. 


To check the validity of SALMON, sample results were compared with those 
obtained using a Monte Carlo model Incorporating the statistical estimation 
methodology of Plass end Kattawar (ref. 9) (PKM). If a photon Is traveling 
directly toward a given detector (counting bln) after any Interaction with 
the water medium, PKM estimates and records the flux expected to reach that 
detector as a function of the optical distance along the current photon path 
to the detector. Sample cases were constructed In which photons at an Incident 
wavelength X - 450 nm were Injected at normal Incidence at an Identical 
spatial location (zero beam spread) Into an Infinitely deep homogeneous water 
medium whose sole constituent was 1.7 x 10^ cells/ml of a common blue-green 
algae ( Anacystls marlnus ). Atmospheric effects were Ignored and the air/water 
Interface was assumed to be perfectly flat. Both SALMON and PKM were executed 
In a dual mode In which photons can be detected not only via single and multiple 
scattering by the algal cells at the Incident wavelength, but also via absorption 
by the chloroplasts contained In the cells and subsequent fluorescent emission, 
with a reduced energy level, at X * 685 nm, the primary fluorescent wavelength 
of chlorophyll a,. These fluoresced photons at X ■ 685 nm are then susceptible 
to scattering and absorption by other algal cells. Scattering and absorption 
coefficients for the algal medium were measured using the cadre of Instruments 
described In Whitlock et al. (ref. 10) and are listed In Table 1 for the two 
wavelengths under Investigation. Absorption coefficients for pure water were 
obtained from Smith and Tyler (ref. 11) while the fluorescence efficiency n 
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• (tht ratio of photons fluoresced to photons absorbad by an algal call) 

* » 450 nm was estimated by combining tha mtasurad absorption coafflclant 

Mlth data provided by Farmer (raf. 12). (Sacond-ordar fluort scene* was 

assumad nagllglbla; hanca n for X ■ 685 nm was sat aqual to zaro.) Tha 
tast medium was purposely chosen to be optically thick to maximize tha 

number of photon Interactions within tha medium* and thus to provjde a 
reasonable upper limit for the computational tte required to execute the 
models. However* SALMON 1$ by no means restricted to optically thick media. 

For the test cases* It was assumed that scattering at both wavelengths* 

and fluorescent emission at X • 685 nm* were Isotropic. For the PKM model* 

the estimated flux exiting the medium within a radius of 1 meter from the 

entry point of the Incident photons was collected In 10 polar-angle counting 

bins (of equal solid angle) ranging up to a maximum emergent polar angle of 

8.1 degrees from zenith. For SALMON* a hypothetical detector geometry was 

defined as an analog to the Innermost collecting bln of PKM* which recorded 

the flux exiting (within a 1 -meter radius) at polar angles less than 2.56 

degrees from zenith. This Innermost bln was found to be the smallest practical 

♦ 

"detector" fleld-of-vlew which would yield results having an acceptable 
variance at a reasonable computer cost using PKM* but Is still an order of 
magnitude larger thar, the fleld-of-vlew commonly associated with LIDAR system 
detectors . 

Comparative normalized histograms of detected signals at A 1 685 nm for 
the test cases using SALMON and PKM are shown In figure 2. The SALMON histo- 
gram shows the distribution of sample means of 20 Independent sets of 1000 
photon trajectories. The PKM histogram stemmed from 36 Independent sets of 
10* photon trajectories; however* since Isotropic scattering and emission were 
specified In the test cases. It was found that the exiting radiance distribution 
was nearly uniform with polar angle* and It was thus assumed that each of the 
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10 counting bins ctuld bo treatod as an Independant collector. Thus, the 
PKM histogram shows the distribution of 360 sample mean detected signals, 
rather than just the 36 collected In the Innermost bln. Even though the 
overall mean estimated signal Is nearly the same f?r the models, figure 2 
shows that the variance In the estimate of the overall mean made using SALMON 

Is only a small fraction of that experienced using PKM. In addition, 

* 

construction of the 36 x 10* photon trajectories using PKM required 15 times 

the computer resources (executing In FORTRAN Extended Version 4 on a CDC 

CYBER 175 computer at the Langley Research Center) expended In construction 

of the 20000 SALMON trajectories. Similar results were obtained for X ■ 450 nm, 

but are not shown In the Interest of brevity. 

To further Illustrate the utility of SALMON, It was executed again for 

the same optical medium, with the exception of substitution of the spectral 

average of the sharply anisotropic measured scattering phase functions for 

Anaeystls marlnus (fig. 3) to govern scattering events at both X ■ 450 nm and 

X ■ 685 nm. In this case, the remote detector was defined as having a horizon- 

2 

tally oriented circular aperture of area .05 cm at an altitude of 150 m, with 
the field of view limited to 6.67 mllliradlans from the vertical (resulting In 
a 1-m radius footprint at the surface). This configuration Is a realistic 
model for helicopter-borne LIDAR remote-sensing syrtem geometries. SALMON 
was executed In Extended Basic on the Hampton Institute POP 11/34 computer 
system In sets of 1000 photon trajectories and the cumulative average detector 
signal (total signal/total number of trajectories) was computed at the conclusion 
of each set. Figure 4 shows a history of the convergence of the cumulative 
average signal for both wavelengths. Fluctuation In the cumulative average 
signal for each wavelength is within +3 percent after only 30,000 trajectories, 
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which required approximately 7.5 hours for execution on the POP 11/34 and 
approximately 1 minute for execution on the CDC Cyber 175. Preliminary 
results obtained using SALMON with much lower concentrations of Anacystls 
marlnus (1.7xl0 5 to 1.7x10* cells/ml) show convergence histories quite 
similar to those In figure 4. Execution tint required for these optically 
thinner cases was somewhat less than that required for the optically thick 
case due to fewer photon Interactions within the Interaction volume. 

This letter has presented a semi-analytic Monte Carlo simulation 
methodology (SALMON), which Is particularly well suited for addressing 
fundamental radiative transfer problems In oceanographic LIDAR, and gives 
a framework for Investigating the effects of a host of environmental factors 
on LIOAR system performance. For the test cases shown, which Include both 
Isotropic and anisotropic scattering and Isotropic fluorescent emission, 
SALMON provides minimal variance, rapidly convergent estimates of detected 
signals for a physically realistic remote sensor geometry. Preliminary 
results for optically thinner media Indicate that performance similar to 
that shown for optically thick media can be expected with SALMON, at a some- 
what reduced level of computational time. 

The authors wish to thank Dr. Wayne Esalas for many hours of fruitful 
discussion, and Mr. David Chrlstllf for Invaluable support In computer 
program translation. We also wish to acknowledge Dr. George Kattawar for 
Inspiring the original thoughts leading up to the development of SALMON. 
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FIGURE CAPTIONS 


Figure 1. 
Figure 2. 

Figure 3. 

Figure 4. 


Remote- sensing geometric model. 

Normalized histograms of sample mean detertor 
signal using SALMON ard PKM. 

Spectrally averaged scattering phase function 
for b^e-green algae. 

Convergence history of cumulative average 
detector signal using SALMON. 
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TABLE 1 - Optical Properties of Test Medium 


Total beam attenuation coefficient, 
m“* 

Scattering coefficient for algae. 

Absorption coefficient for algae. 

Absorption coefficient for water, 
m"* 

Fluorescence efficiency 


A ■ 450 nm 

A ■ 685 nm 

16.01 

8.20 

13.70 

7.00 

2.28 

0.75 

0.03 

0.45 

0.00254 

0.0 




Normalized frequency of occurrence 
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Sample mean detector signal, relative units 

Figure 2. - Normalized histograms of sample mean detector 
signal using SALMON and PKM. 
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Cumulative average detector signal, relative units 
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(b) X = 685 nm 

Figure 4. - Convergence history of cumulative average 
detector signal using SALMON. 
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Chapter 2 


Sensitivity of Airborne Fluorosensor Measurements 
to Linear Vertical Gradients in Chlorophyll Concentrations 


D. D. Venable, A. R. Punjabi and L. R. Poole 



Introduce Ion 


I 2 

We employ a semianalytic Monte Carlo radiative transfer model (SALMON) ’ 
to investigate the effects of vertical inhomogeneities in phytoplankton con- 
centration inTlear sea water on airborne laser-induced signals due to 
chlorophyll fluorescence and water Raman scattering. In this model, the de- 
tector signal is evaluated through a combination of stochastic and analytic 
techniques making it possible to model lidar systems with realistic geometric 
constraints and at the same time reduce required computer resources. 

A Layered approach is used to characterize changes in phytoplankton 
concentration, and thus optical parameters, as a function of depth. This 
approach allows media to be represented by combinations of homogeneous layers 
of various thicknesses. We assume that appropriate optical parameters can be 
defined for, and that horizontal homogeneity holds within each layer. 

The geometry of the physical set-up used in the model is shown in Fig. 1. 

1 2 

The geometric model is the same used in previous studies * except that the 
medium is divided into 21 layers, where layer l is always the atmosphere. 

The clear water medium, containing chlorophyll is divided into 20 layers. 

The concentration of chlorophyll in each layer must remain constant, but con- 
centrations may vary from layer to Layer. Our current modeling effort considers 
only linear variations in concentration and assumes that the airborne laser 
platform is at an altitude sufficiently Low that atmospheric effects can be 
neglected. 

In Fig. I, q is the termination point of the j** 1 trajectory of the photon 
packet, € is the direction of propogation of the packet, and 0’ is the angle 


through which photons must be scattered at point q in order to be collected 


by the detector within the solid angle AO. A is the area of the detector 

ap 

aperature and the detector is Located at the laser source, 150m abov« the 
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water surface# The detectors field of view and the interaction volume are 

represented by broken lines in the figure. The distances dj, d^, • • # refer 

to the distances that a photon packet must travel through layer 2, layer 3, 

. . . to reach the surface, assuming it had scattered through an angle at 

point q. The refractive index of the water medium is assumed to be the same 

for each layer. Thus, no refractions occur at the boundaries within the medium. 

Each layer thickness may be set to any desired value. Thicknesses were 

selected so that for a given surface concentration, the depth of the bottom 

of layer 20 is approximately two to three times the maximum signal-integrated 

depth, max * obtained from homogeneous results for the three wavelengths 

involved. The values we used for Z n _ were taken from Ref. 3 for homogeneous 

90, max 

concentrations that were the same as the surface concentrations used in the 

layered model. The signal depth, Zg Q , is that value above which 907. of the 

3 

total signal observed at the detector originates. The thickness of the 21st 
layer is set to a large enough value that photons do not reach the bottom of 
the medium. Optical properties for each layer are based on the layer’s chloro- 
phyll concentration. The concentration for the medium can be denoted by C(z), 
where z is depth in the medium. The downward direction is taken to be positive 
and the surface concentration is given by C Q . The gradient is defined to be, 
g, where 


1 dC(z) 

8 = C dz 
o 

and is given in units of m”*. In the layered model, the concentration 
assigned to the k tl1 layer is the averaged values of C(z) at the top and 
bottom of the layer. Concentration in the 21st layer is assigned so that the 
difference in concentrations for the 21st and 20th layers is the same as the 
difference in concentrations for the 20th and 19th layers. Computations 
for negative gradients of -107. m -157. m”* , and -207. m""* - were not made for 
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surface concentrations less than 1 pg/l, since, in these cases the concen- 
tration of chlorophyll would reach zero at a depth less than Zg Q maj< . These 
concentration profiles would thus deviate from the Linearly varying case and 
could not be compared to ocher cases considered in the simulation. 

Special consideration must be given when the photon packet crosses the 
* 4 

interface between two layers. For example, consider the case where the 
photon packet travels from layer k to layer k t 1. We select the distance 
the photon packet travels, 5^, from an appropriate distribution function.^ 
The value of 6^ is the distance the photon packet would travel in a homo- 
geneous media, that is, if the packet remained in the k C ^ layer. However, 
if a layer boundary is crossed, the additional distance that the photon 
packet travels in layer k ± l is given by 
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k ± l 




( 2 ) 


k k t L 

where is the actual distance traveled in Layer k and and arc 

the total attenuation coefficients in Layers k and k ± l, respectively. 

In this study, the common coastal zone dinof lagel late Prorocentrum 

minimum was assumed to be the only constituent in a medium of clear ocean 

water. Pertinent optical parameters used in the simulations were the same as 

those used in Ref. 3 and are given in Table L. The excitation source was 

assumed to be a laser with an emission wavelength of 480nm. Raman scattering 

at 574nm ond fluorescence at 605nm were assumed to occur at sharp spectral 

lines; i.e., no effects of spectral bands for these- events were considered. 

The scattering phase function for chlorophyll given in Fig. 2 that was used 

3 

for this study is the same as that used by Poole and Esaias. Ref. 3 gives 
appropriate discussions on the origins of the data in Table l and Fig. 2. 
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Results and Discussion 

We have focused our study on the water Raman normalized Fluorescence 
signal, R = H p /H R , where H p is Che signal at 685nm received at the detector 
and H R is the signal received at 574nm. In particular, we investigated R as 
a function of vertical concentration gradients for surface chlorophyll concert- 
trations in the range from 0.01 to 20 Mg/ 1 * Simulation results for gradients 
of ± 207., ± 15%, ± 10%, and ± 5% m” 1 are shown in Fig. 3 where R is plotted 
as a function of surface concentration for each pair of gradients. The cor- 
responding plot for a homogeneous medium with a concentration equal to the 
surface concentration of the analogue inhomogeneous media is shown on each 
graph for comparison. Error bars are within the size of the symbols used for 
plotting. The functional relation between R for a particular value of g, R , 

O 

and surface concentration, C q , for the data shown in Fig. 3 is of the form 

R (C ) = a C b , (3) 

goo* 

where a and b are constants for a given gradient. It she uld be observed that 
since the values of b are near unity (ranging from l.l - 1.3), the data arc 
well approximated by a linear relationship as indicated in Ref. 3. This is 
particularly the case when only the larger values of surface concentration 
are considered. 

Absolute calibration of airborne fluorosensor measurements is practically 
impossible due to factors such as variation in phytoplankton fluorescence 
efficiency (or cross-section), both among the various species and also with 
environmental conditions' (ambient light, nutrient levels, etc.). As a result, 
measurements are calibrated empirically by comparison with iri situ sampling 
at or near the water surface; obviously, vertical inhomogeneities in chlorophyll 
concentration could lead to bias in the empirical calibration. To establish 
ranges of surface concentrations and/or vertical gradients in which calibration 
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by surfaca samples would be Inappropriate, we have applied a one-tall t-teat 
to determine at the 95% confidence level If R for a given vertical gradient 
Is significantly different chan the value of R computed for the analogous 
homogeneous medium. We assume that the results for the homogeneous and In- 
homogeneous cases represent Independent samples. Approximately ten runs are 
made to obtain values for R for each gradient and for each surface concentra- 
tion. The computer simulations are terminated when the variances In Hp and 

H_ fall below 5% of the values of H- and H D and at leajt. 30 thousand photon 
R r r 

packets have been sampled. We assume that our two secs of samples have 
common variances for the purpose of the t-t^st. Results of the t-test are 
given in Table 2. For media with high surface concentrations of chlorophyll 
(> 10 pg/l), the fluorosensor measurements are Insensitive to vertical con- 
centration gradients. For cases in which the surface concentrations are 
less than l vg/l, fluorosensor measurements were found to be sensitive to 
even the mildest gradients tested. For the intermediate range of surface 
concentrations ( l - 10 jig/l), the sensitivity of fluorosensor measurements 
to gradients was found to generally increase as the magnitude of the gradient 
increased. Significant differences (at the 95% confidence level) between 
the mean values of R for corresponding homogeneous and inhomogeneous cases 
occurred at surface concentrations of 1, 6, 6, and 10 pg/l for gradients of 
t 5%, ± 10%, ± 15% and ± 20% m“ l , respectively. 

Concluding Remarks 

We have extended a semianalytic Monte Carlo radiative transfer simula- 
tion model for airborne laser f luorosensors to allow for Inhomogeneities in 
the vertical distribution of consitutcnts in clear sea water. This article 
has presented results of the simulations for lineary varying step concentra- 
tions of chlorophyll. The results indicate that statistically significant 
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differences can be seen, under certain conditions, in the water Raman normal- 
ised fluorescence signals between nonhomogeneous and homogeneous cases. 
Generally, remote sensor users must exercise care when making empirical 
calibrations of airborne f luorosensors. Tha SALMON model has been used to 
show that fluorosensor measurements are sensitive to certain ranges of sur- 
face concentrations when the chlorophyll concentrations vary lineary with 
depth in the water medium. 
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Figure Captions 


Figure L - Geometric arrangement of the f luorosensor system modeled by 
SALMON. 

Figure 2 - Spectrally averaged scattering phase function for Prorocentrura 
minimum . 

Figure 3 - Graphical representation of the Raman normalized fluorescence 
signal, R = Hp/Hg, vs surface concentration for gradients of 
(A) ± 207. m"\ (B) ± 157. nf l , (C) ± 107. nf l , and (D) ± 57. m" L . 

A - positive gradient, V - negative gradient, and 0- homogeneous 


distribution. 



Lise of Tables 


Table l 
Table 2 


- Ope lea l Parameters 

- Confidence Limies for Deeeceing Significant Differences Between 
Mean Values of R for Corresponding Homogeneous and Inhomogeneous 
Cases 
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Table 1; Optical Parameters 



480 

Excitation 

Wavelength, 

574 

Raman 

nm 

685 

Fluorescence 

Absorption coefficient for P. minimum* 
m” V (Mg/liter) 

0.0582 

0.0382 

0.036 

Scattering coefficient for P. minimum, 
m" l /(yg/liter) 

0.291 

0.231 

0.202 

Fluorescence efficiency for P. 
minimum 3 

0.028 

— 

— 

Absorption coefficient for clear ocean 
water, m“^ 

0.0176 

0.0911 

0.475 

Scattering coefficient for clear ocean 
water, m“^ 

0.0034 

0.0017 

0.0007 

Raman scattering coefficient, m”* 

0.0018 

0.0009 

0.0004 


a Fluorescence at 685 nm considered only for excitation photons at X = 480 nm. 
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Table 2s Confidence Limits for Detecting Significant Differences Betvtten 
Mean Values of R for Corresponding Homogeneous and Inhomogeneous Cases 


Surface 

Gradient 

Gradient 

Gradient 

Gradient 

Concentration 

ng/i 

± 207. nf 1 

± 157. nf 1 

± 10% nf 1 

± 5% nf 1 

0,01 

99 s- 

99+ 

99+ 

99+ 

0.03 

99+ 

99+ 

99+ 

99+ 

0.1 

99+ 

99+ 

99+ 

99+ 

0.3 

99+ 

99+ 

99+ 

99+ 

1 

99+ 

99+ 

99+ 

98 

3 

99+ 

99+ 

96 

89 

6 

98 

99 

95 

90 

10 

95 

92 

88 

54 

15 

86 

51 

78 

58 

20 

79 

77 

67 

58 
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Chapter 3 


SALMON Modifications for Non-hornogeneous 
Media 


A. R. Punjabi 



Introduction 


We have employed the semianalytic Monte Carlo radiative transfer model 
(SALMON)* to investigate the effects of inhomogeneities in the phytoplankton 
concentrations in the clear sea waters on the lidar signals. In this model, 
the detector signal is evaluated through a combination of stochastic and analytic 
techniques. This makes it possible to model lidar systems with realistic geometric 
constraints such as small, finite detector fields of view. Inclusion of analytic 
estimates for various stochastic processes reduces the computer resources re- 
quired to use the model. Convergence in the SALMON model is virtually independent 

of collector size. Significant statistical advantages are gained by splitting 

2-4 

true Raman photons and modifying virtual interactions. 

The geometry of the physical set up used in our investigation is shown in 

2 

Fig. 1. A ap is the area of detector aperture (0.071 cm ). Layer l is air. Broken 
lines in the layer l represent the detector field of view (0.0025 radians). 

However, the SALMON simulation code can handle l to 5 fields of view. The detector 
is located at a height of 150 m above the top of water surface which is perfectly 
flat. The clear water medium containing phytoplankton is divided into layer 2 
through layer 21. Concentration of chlorophyll is constant in each layer but 
varies linear.y from layer to layer as one advances from top of layer 2 downwards. 
For a given surface concentration layers 2 to 20 have equal thicknesses. The 
bottom of layer 21 is 999 ms below the air-water interface for all values of sur- 
face concentration. The thickness of 2lst layer is set to such a large value to 
ensure that the photons do not reach the bottom of medium. Thicknesses of layers 
2 to 20 is large (small) for small (large) surface concentrations. Each layer is 
completely homogeneous in the chlorophyll concentration. No lateral variations 
are allowed in chlorophyll concentration. The linearly varying profile of chloro- 
phyll concentration in water medium is approximated by a step function profile. 
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In Fig. I, q is the terminal point of J th trajectory of the photon packet, 
q is its direction of travel. 0' is the angle through which the photon at point 
q must be scattered in order to be collected by the detector within the solid 
angle Afl. Broken lines in layers 2 to 21 define the boundaries of the interaction 
volume . ^ 

Fig. 2 is the flow chart for the SALMON simulation code. The number of events 
N q per sample is limited to 1000, and number of samples is restricted to a maximum 
of 120 with a minimum size of 30 samples. 


Input Parameters and Other Pertinent Quantities 

Optical parameters used in our computations are shown in Table l. In this 

study, the common coastal zone dinof lagel late Prorocentrum minimum is assumed to 

be the only constituant in a medium of clear sea water. Pertinent optical parameters 

used in the simulation are the same as those used in Ref. 2. The spectrally averaged 

volume scattering phase function for Prorocentrum minimum, which we have used in 

2 

this study, is the one used by Poole and Esaias , and is shown is Fig. 3. Ref. 2 
gives appropriate discussions on the origin of data in Table l and Fig. 3. 

The surface concentration of chlorophyll is denoted by C q . The concentration 
for the medium is denoted by C(z), where z is the depth in the medium. The down- 
ward direction taken as positive. The gradient in chlorophyll concentration, 
denoted by g, is defined to be 


1 dC(z) 
8 “ C dz 


( 1 ) 


o 

and is given in units of m - ^. In the stratified model, the concentration assigned 

A* L 

to k layer is the averaged value of C(z) at the top and bottom of the layer. 
Concentration in the 21st layer is assigned so that the difference in concentrations 
for the 2ist and 20th layer is the same as the difference in concentrations for the 
20th and 19th layers. 
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Z 9o is defined as the depth from which 907. of detector signal at a particular 
wavelength emanates for a homogeneous concentration corresponding to C Q . ^ 9 Q >max 
is defined as the maximum of the for laser wavelength ^ )» z 9 ^ for the 

Raman wavelength (Z 9Q ^ ) and Z 9Q for the fluorescent wavelength (Z^ ^ i * e ’ 

Z 90,max ( V C o ) * m * ( Z 90, X L Z 90, X R Z 90, X p (C o ) } (2) 

Note that Z Art depends on the laser excitation wavelength X. and the value of 
vo,max *-< 

surface concentration C Q . Z 9 ^' s used in equ. (2) arc for the homogeneous medium 

and are taken from Ref. 2. Fig. 4 shows vs. for X^»532 and 480 nm. The 

excitation source is assumed to be a laser with emission wavelength of 480 nm or 

532 nm. Raman scattering at 574 nm and fluorescence at 685 nm are assumed to occur 

at sharp spectral lines; no effects of spectral bands for these events are considered. 

For a given surface concentration, the depth z q in the water medium below 

which the concentration C(z) vanishes for negative gradients g, depends upon g 

and C , Z is defined as 
o o 


C(z > z ) £ 0 , g < 0 . (3) 

o w 

z may be smaller than Z nn for a given „xcitation wavelength when g is a 
o yo,max 

sufficiently large negative number. Table 2 shows Az and Z QA for various 

7U,n la A 

values of surface concentrations G q where Az is the thickness of layers 2 through 

20. Table 2 is for excitation wavelength x^=480 nm. This table also indicates 

the depth z q for values of gradients g = -207. m -157. m *, 107. m * and -57. m 

Solid demarcating line in the last 4 columns of this table separate the C Q regime 

for which z < Z n „ from the C regime for which z > Z n . . For example, 
o 90, max o & o 90, max 

for a gradient of -207. m~ and surface concentration below 1.0 pg/l, z q is smaller 

than Z rt _ . C* denotes the value of C for which z is smaller than Z n _ 

90, max o o yo„max 

for a given gradient for all C q <_ C*. Table 3 shows C* for various values of 
gradients for *^=480 nm. 
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In this study, wo have considered the range of surface concentrations from 

0.01 pg/l to 20 pg/l, while the range of gradients is -207. nT* to 20% rrf l . Table 

3 is based on these ranges of values of C Q and g. It must be pointed out that 

for excitation wavelength A L »532 nm, U 9o>max l is always smaller than JzJ for the 

ranges of C o and g considered in our study. This distinction between the two 

excitation wavelengths X L »480 nm and 532 nm is of considerable significance in 

this study. For A£«480 nm, we can see from Fig. 4 that Z^q ^ is large compared 

R 

to Z n „ , or Z nn . for small values of surface concentration. Thus, Z Q _ * 

90, Xp 90, A^ yu,max 

Zq_ for X,®480 nm. Thus, in this case, when surface concentrations are equal 

“v» Xp L( 

to or less than C*, the detector receives 907. or more of signal at fluorescent 
wavelength while less than 907. of signal is received at Raman wavelength. This 
asymmetry between signals received at the Raman and fluorescent wavelengths becomes 
manifested when we Txamine the Raman normalized fluorescence signals as function 
of surface concentration for a given gradient when C q > C*. We will observe that 
such points of asymmetry do not give good fits to the simulation results. While 
we will also observe that no such difficulty arises for excitation wavelength 
A L »532 nm. 

In the following discussion, we denote the detector signal received at the 
Raman wavelength by H^, while the detector signal received at the fluorescence 
wavelength is denoted by Up. Both Hp and are functions of the excitation wave- 
length A l , the gradient in chlorophyll concentration g and the surface concentration 


C q . The Raman normalized fluorescence signal is denoted by R, where 

RU l , C o , g) = H f (X l , C o , g)/H (X l , C 0 , g) 


(4) 


For a fixed surface concentration and gradient g, ten computer runs are 
made. Outcome of a particular run is accepted if 957. of the samples have 107. or 
less deviation in signals Hp and H^. R is then calculated for each successful 
computer run from eqn. (4) and the standard deviation in R, a is calculated 


36 



from the relation 


°R ' R «*S,/»J> * <°H 

r K 

where and o R are standard deviations in H p and H R rospeetivoLy 

F R 

values of R and standard deviation a R are the weighted averages of 
given by 


<R> = 


n 

fi (Rj / 
i=i 



n 

B (1/ 

i=i 


R 


2 

i 


and 


n 

<<J >ss ± {1/ E 

K i*i 


( 1 !«l >>* 


whore 

n < 10 

For the remainder of this report, R will represent <R> and a R 

< °R > - 


(>) 

. Then the mean 
P and o R , 

( 6 ) 

( 7 ) 

will represent 
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Results and Discussion 


Tables 4 and 5 show R and o R values as functions of surface concentration 
C and gradient g for excitation wavelengths X,*480nm and 532nm respectively. 

Let us first consider the simulation results for A L *480nm. We first attempted 
to fit the simulation results of R vs C o , for a given g, by a straight line with 
zero intercept. However, this fit was ouite poor for low surface concentrations. 

The reason for requiring a zero intercept is that when chlorophyll concentration 
vanishes, the fluorescence signal NL will vanish. A straight line with zero 
intercept can bo represented by 

R(C q , g) s* m(g) C q (8) 

where m is the slope and m varies with gradient. o m (g) ^ ihc standard deviation 
in m. It must be mentioned that for A,~480nm, we fit; only those points for which 

IJ 

z > 7 nn . If we attempt to fit all points (including those for which 
o 90, max r 

z < Z nn ), the overall quality of fit deteriorates. When we plot log R vs 
o 90, max ’ ^ 

log for a given gradient g, a straight line gives a good fit. Figs. 5(a) - (d) 
show this power law fit for A L ~480nm. The power law relation is 

R(C o , S) = A C® (9) 

where the scaling factor A and the exponent B are functions of g alone for a given 

A, . cr. and o_ represent standard deviations in A and B respectively and are 
LA B 

functions of g. 

Figs. 6(a) - (d) show straight line fits (with zero intercept) for the A^=532nm 
simulation results. Figs. 7(a) - (d) show power law fits for the same excitation 
wavelength. It is clear that a straight line fit with zero intercept is more 
satisfactory in case of green laser than in case of blue laser. However, in both 
cases, the power law fit is superior to the straight line fit. 
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Tables 6(a) and (b) show the fitting parameters m, A & B and standard 
deviations o ^ and as functions of gradient for and 532nm respectively. 

Note that both A and m decrease with gradient while B increases as the gradient 
decreases for both excitation wavelengths. Uncertainties o^ and for both 
excitation wavelengths become larger as the gradient becomes smaller for positive 
vaLues through negative values. Fitting parameter B for the Xp»532nm case is 
very close to unity for all values of g. Also, m and A are quite comparable for 
Xp=532nm. The standard deviations, and for Xp=532nm are smaller than those 
for \ps480nm for ail values of g. 

2 2 

Now we define two goodness-of-f iu parameters, X^ and X 2 , as 

x l • R £ (R fu - R i>'»i 2 * < l0 > 

and 



N 

Z (R 
i=i 


fit 



(ID 


where is given by eqn. (8) for a straight line fit with zero intercept and 

2 

by eqn. (9) for a power law fit. weighs each equally. In eqns. (10) and 

2 

(11) above, N is the number of points. Thus X^ weighs R for small and large 

2 

concentrations equally. While X 2 weighs R for large heaviLy compared to R 

2 

for small C () , for a given g. We also define as ratio of X 2 for straight line 

2 2 
fit to X 2 for power law fit, and P^ as the ratio of X^ for straight line fit to 

2 

X^ for power law fit, i.c. 


and 


P 


2 


2 

X 2 » straight line fit' 
~2 

X 2 , power law fit 


( 12 ) 


VI 
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X, , straight line fit 

'l = "i ; 7 

power law fit 


( 13 ) 


Tables 7(a) and (b) give Xp X 2 * P^ and P 2 for both straight line fit and 
power law fit for X,»480 and 532nm, respectively. 

From Tables 7(a) and (b) we see that, for high surface concentrations, a 
straight line fit is as good as a power law fit and it is all the more so in the 
case of the green laser. While if we consider ail surface concentrations, a 
power law fit is definitely far superior to a straight line fit for both blue and 

V 

green lasers. Thus we can definitely conclude that a power law relationship given 

by eqn. (9) is the better - alytic representation of the computational results. 

2 2 

Let us examine the average values of and for both straight line fits 

and power law fits for X^=480 and 532nm as 

2 12 
<X 1 > = 9 g X l,g 


2 12 
<X 2 > = 9 l X 2,g 


(14) 


where the summation is over various vaLues of gradient g. Table 8 shows these 
average values. This table clearly shows that a power law fit is better and that 
computational results for green laser fit better. 

If we fix the surface concentration and average 0 ^ over all values of 
gradient g, we get an average value for c_, denoted by <0 O > , as a function of 

K Kg 

C q . This can be expressed as 


< a_ > = 7; Z o_ 

R S 9 g R >g 


(15) 


where again the summation is over various values of gradient g. Table 9 lists 
<0p> for both excitation waveiengLhs. We observe from this table that errors in 
R, the Raman normalized fluorescence signal, are smaller at low values of surface 
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concentration and also that for C 0 »C'.Cl to 3 Mg/ 1, <0 R > g art ‘ smaller for green 
laser than for blue laser and for C q > 3 pg/1, the opposite is true. 

In Figs, 8(a) and (b) we show power law curves and computational points for 
gradients g » 20% m~* # 07. m"* anJ -207. m - * for bLue and green lasers respectively. 
Figs. 9(a) and (b) show power law fitting parameters A and B vs gradient g for 
blue and g ion lasers respectively. 

Fig. 0 and Tables 4 and 5 show that signal ratio P. is higher (lower) for 
positive (negative) gradients as compared to the homogeneous case for all values 
of surface concentrations considered in this study. The exponent in the power 
law fitting function decreases as gradient g increases from -207. m" * through 
0% nf l to -207. nf l . 

The upper limits on the range of surface concentration values beyond which 
error bars on R(C t) , g) and RfG^, g=o) and/or R(C q , g=o) and R(C q * «g) overlap 
are denoted by and are shown in Table 10. This table shows that blue and green 
laser results are similar so far as overlapping of signals from homogeneous media 
with the one from inhomogeneous media is concerned. As is expected, the upper 
bound on the range of surface concentrations within which signals do not overlap 
reduces with decreasing magnitude of gradient. 

We have also performed one tailed T tests with unequal variances to determine 

distinguishabil ity of R(C o> ± g) from R C , g=0). Tables 11(a) and (b) show the 

results of these tests. The tables show the probabilities a + or a_ that signal 

R(C o » - g) is greater (or less) than R(C q , g=o). a is the smaller of a + and a 

for a given C q . We see that for lower surface concentrations, signal R for a 

given positive or negative gradient is distinctly higher or lower than R for zero 

gradient. Thus, for smaller surface concentrations, R from inhomogeneous media can 

easily be distinguished from the R from a homogeneous media with the same surface 

concentration C . 

o 


We define C^ s as Che upper limit: on beyond which sigral. R from an 
inhomogeneous media with a given magnitude of gradient, |g|, can not be distin- 
guished from R from a media with g=o, with a > 95%. Table 12 lists C^ s as a 
function of |g| . 

From Tables ll and 12, we see that the range of surface concentrations 
for which signals can bo, resolved between homogeneous and inhomogeneous media 
becomes smaller as the magnitude of gradient becomes smaller, as can be expected. 
C D for green laser are better than C _ for blue laser. Thus green laser will 
resolve signals (between homogeneous and inhomogeneous media) better, i.e. for 
wider range of surface concentrations. 
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Conclusions 


We have employed the semianalytic Monte Carlo simulation methodology towards the 
following ends: 

(i) To determine how the lidar signals from an inhomogeneous media 
differ from ones from homogeneous media. 

(ii) To determine an analytic representation of simulation results. 

(iii) To compare the performance of blue and green lasers. 

In the SALMON simulations we performed, the linearly varying vertical 
chlorophyll concentration profiles have been approximated by step function profiles. 
For this reason, we refer to the simulation model as the stratified SALMON model. 

The two parameters which we use to specify the inhomogeneous media are C q , the 
surface concentration and g, the gradient in concentration. 

The simulations have shown that the Raman normalized fluorescence signal, R, 
is generally higher (lower) for higher (lower) surface concentrations C q and 
gradients g. R is high for positive gradients compared to zero gradient. Opposite 
results are true for negative gradients in chlorophyll concentrations. It is 
easier to dintinguish media having different gradients at lower values of surface 
concentrations. At high values of surface concentrations, the differences in the 
signals frommedia having different gradients overlap and can not be distinguished. 
The range of surface concentrations and the upper bound of this range for which 
signals from media having a fixed non-zero magnitude of gradient can be differen- 
tiated from those which emanate from homogeneous media reduces as the gradient 
becomes smaller. In short, lower concentrations and high gradients yield better 
resolution and distinguishability of signals. 

The normalized signal R and the surface concentration C q can be approximated 
by an analytic expression of the form 

R <c o , g> = A (g) c o B <&> 
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( 16 ) 



In this power law relationship, the scaling factor A and the exponent B are 
functions of g alone and both show a definite trend in their relationship to g 
(see Fig. 9). This power law fit is better for small C q *s as compared to high 
C o 's. Since error bars on A and B do not overlap beyond more than 10% m ^ 
difference in gradients, this analytic representation is good for media having 
gradients differing by more than 10% m“^. 

The relative ti in- in R given by cr R / R , where o R is the standard deviation in 
R, are smaller for green laser excitation than those for blue laser excitation. 
Power law fits for A L =532nm are superior to those for A L ~480nm. Green laser 
e;<citation can resolve greater range of surface concentrations than blue laser 
excitation. Thus the computational results show that results for green laser are 
more satisfactorily standarizable (by an analytic expression) and more widely 
applicable (for the purpose of distinguishabi lity of signals) than those for blue 
laser. 
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TABLE 1: Opt I rzl PoroMtltrt 


Wavelength , nis 


480 574 685 

Excitation Ron an Fluorescence 


Absorption coefficient for P. minimum, 
/m/Cpg/ 1 ) 

Scattering coefficient for P. minium, 
/m/Cjjg/l) 

0.0582 

0.0382 

0.038 

0.291 

0.231 

C .202 

Fluorescence efficiency for P. 
a 

minimum 

0.028 



Absorption coefficient for clear ocean 
water, /m 

0.0176 

0.0911 

0.475 

Scattering coefficient for clear ocean 
water , /m 

0.0034 

0.0017 

0.0007 

Raman scatternlg ‘coeff lelent,/m 

0.0018 

0.0009 

0.0004 


a 

Fluorescence at 685 ran considered only for axel tot Ion wavelength at 480 run. 
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TABLE 2 


EXCITATION WAVELENGTH-480 dm. 


c 

0 

C^o/I) 

A 

Cm) 

Z 

90, a ax 
Cm) 


d 

0 

Cm) 



- 20X 

9 

C/m) 

- 15 % 

s 

1 “ 

- 5X 

.01 

2.0 

16.0 

6.0 

8.0 

12.0 

22.0 

.03 

2.0 

17.0 

6.0 

8.0 

12.0 

23.0 

.1 

2.0 

13.5 

6.0 

8.0 

12.0 

22.0 

.3 

2.0 

9.0 

6.0 

8.0 

fn.e 

22.0 

1. 

1.0 

5.0 

xr 

8.0 

11.0 


3 . 

1.0 

3.0 





6.0 

0.1 

2.0 





10.0 

0.1 

* 1.8 





15.0 

0.1 

1.6 





20.0 

0.1 

1.4 
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TABLE 3 


EXCITATION UAVELENGTH-488 no. 

» 


c 

9 

Cjjq/\) 

C/m) 

B.3 

-20% 

0.1 

-10X 
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Exc I lotion wavelength St. I Ine fit Power low fit 
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0.0769 0.0448 0.0020 0.0390 


56 


ORIGINAL F:\w- t3T 
OF POOR QUALITY 


TABLE 9 


Surface Concentrat Ion Average Error <S > 

R r 


C Excitation Wave length 

0 


C/jq/\> 489 nn S32 m> 


8. 91 

8.8118 

8.8898 

9.93 

9.9116 

8.8892 

8.1 

8.8125 

8.9194 

9.3 

8.8177 

9.8152 

1.9 

8.9389 

8.9254 

3.8 

9.8472 

9.8453 

6.8 

8.8587 

8.8531 

* 18.9 

8.9473 

9.8548 

IS. 9 

8.9541 

8.8621 

28.9 

8.8673 

9.9765 


57 



C/iO 


20X 

15X 

10X 

5X 


ORIGINAL PAGE &S 
OF POOR QUALITY 


TABLE t0 


C 

OL 

Cfig/I > 

Excitation wovalongth 


480 nn 

532 nn 

10.0 

10.0 

10.0 

10.0 

6.0 

1.0 

1.0 

1.0 


58 




BOlfiflltt- 


ORIGINAL PAGE IS 
OF POOR QUALITY 


TABLE 1 1 Co) 


C Gradient 

0 

WD C/iO 

20X -20X 15X -15X I0X -I0X SX -5X 


00.01 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

00.03 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

3 

A 

00. t 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

00.3 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

1.0 . 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

98 

>99 

98 

3.0 

99 

>99 

99 

>99 

>99 

>99 

90 

>99 

90 

89 

98 

89 

0.0 

98 

>99 

90 

99 

>99 

99 

95 

>99 

95 

90 

>99 

90 

10.0 

95 

90 

95 

92 

99 

92 

98 

88 

88 

81 

54 

54 

15.0 

98 

80 

88 

87 

51 

51 

80 

78 

78 

58 

68 

58 

20.0 

79 

88 

79 

77 

89 

77 

67 

82 

87 

58 

68 

58 



CD<n oou- 




W- 


ORIGINAL PASS m 
OF POOR QUALITY 


table: i tcb) 


C Gradient 

0 

C»iq/1) C /■»> 


26X -20X 15X -15X 10X -10X 5X -5X 



*** 


< 



4 



«c 



•C 

1.01 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

1.03 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

1.1 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

1.3 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

>99 

.0 

>99 

>99 

>99 

>99 

>99 

^99 

>99 

>99 

>99 

99 

>99 

99 

1.0 * 

98 

>99 

98 

98 

>99 

98 

96 

93 

>99 

89 

90 

89 

1.0 

>99 

>99 

>99 

98 

>99 

98 

98 

90 

90 

88 

71 

71 

1.0 

98 

>99 

98 

99 

>99 

99 

96 

98 

98 

92 

98 

92 

.0 

92 

65 

65 

79 

68 

68 

85 

79 

79 

53 

88 

53 

1.0 

81 

79 

79 

77 

80 

77 

68 

82 

68 

69 

76 

69 


60 


ORIGINAL Pt/Al i‘J 
OF POOR QUALITY 


TABLE 12 
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Figure Captions 


Uhl-K •• - ' 7 
OF POOR Q--'- 


.Is * 


1. Geometry of the physical set up used in our investigation. 

2. Flow chart for SALMON simulation code. 

3. Spectrally averaged volume scattering phase function for Prorocentrum minimum. 

4. Z^q vs. C q for X^=480 and 532 nm. 

5. Power-law fits for X^-480 nm. 

(a) g = 20% m ~ l , 0% nf 1 and -20% m“ L , 

(b) g - 15% m" 1 , 0% m" 1 and -15% m“ L , 

(c) g = 10% m”^, 0% m - ^ and -10% m S 

(d) g = 5% m”* - , 0% m”^ and -5% m 

6. Straight line (with zero interapt) fits for X L =532 nm. 

(a) g = 20% m~ L , 0% m" 1 and -20% m , 

(b) g « 15% m~ l , 0% nf l and -L5% nf 1 , 

(c) g = 10% nf 1 , 0% nf 1 and -10% m , 

(d) g - 5% m” 0% m~ * and -57„ m 

7. Power-law fits for X,=532 nm. 

-1 L -l -1 

(a) g = 20% m , 0% m , and -20% m 

(b) g = 15% m *■, 0% m ^ and -15% m * , 

( c ) g = 10% nf 1 , 0% nf 1 and -10% nfl, 

(d) g = 5% nf 1 , 0% nf 1 and -5% m. 

8. Power-law fits curve and computational points on a linear plot for g== +20% m 
0% m _1 and -20% m 

(a) for X^=480 nm, 

(b) for X^=532 nm. 

9. Power-law fitting parameters A and B vs. gradient g. 

(a) for X^=480 nm, 

(b) for X^=532 nm. 
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Appendix I : Comparison of Stratified SALMON and Original Homogeneous Models 


The original homogeneous SALMON model is operative on the NASA LaRC Cyber 
computer system. The stratified model is operative on the Hampton Institute VAX 
11/780 computer system. The comparisons here are made for the zero gradient 
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